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We develop a theory of soliton spiraling in a bulk nonlinear medium and reveal a new physical 
mechanism: periodic power exchange via induced coherence, which can lead to stable spiraling and 
the formation of dynamical two-soliton states. Our theory not only explains earlier observations, but 
provides a number of predictions which are also verified experimentally. Finally, we show theoretically 
and experimentally that soliton spiraling can be controled by the degree of mutual initial coherence. 



Self-guided optical beams (or spatial solitons) have 
attracted substantial research interest in the last three 
decades Jlj]. Although interactions between two- 
dimensional (2D) solitons in Kerr and non-Kerr media 
have been studied extensively, only the recent discoveries 
of stable three-dimensional (3D) solitons in different non- 
linear bulk media (|] initiated an experimental study of 
fully 3D interactions of solitary waves. Recently, exper- 
iments demonstrating non-planar interaction and spiral- 
ing of spatial solitons in a photorefractive medium have 
been reported |J. However, in spite of earlier predictions 
of non-planar soliton interactions the experimental 
results || have not been explained theoretically so far. 
Also, it has been shown [|| that coherently interacting 
solitons do not allow any stable spiraling, in sharp con- 
trast with the experimental observations j3|. The funda- 
mental question remains: Is soliton spiraling possible at 
all as a stable dynamical regime of soliton interaction? 

In this Letter we develop, for the first time to 
our knowledge, a theory of soliton spiraling in a 
photorefractive-type nonlinear bulk medium. We derive 
an analytical model describing stable soliton spiraling 
and predict a number of new effects in soliton interac- 
tions, such as an induced coherence and control over 3D 
interactions, which we verify here experimentally, using 
experimental setup similar to that reported earlier [||. 
Importantly, our analytical model and numerical simu- 
lations show that interacting-spiraling solitons conserve 
angular momentum. We believe that this result is a core 
foundation for future research on 3D soliton control, re- 
sembling the conservation of linear momentum in the in- 
teraction of more conventional (l+l)-dimensional soli- 
tons §. 

We consider incoherent beam interaction in an 
isotropic saturable nonlinear medium described by two 
coupled normalized nonlinear Schrodinger (NLS) equa- 
tions 
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where u and w are the beam envelopes, z is the prop- 
agation distance; = d 2 /dx 2 + d 2 /dy 2 accounts for 
the diffraction in the transverse (x, y) plane. This sys- 
tem, in the 2D case (i.e., for = d 2 jdx 2 \ gives rise to 
incoherently-coupled soliton pairs j?J and to incoherent 
collisions || which have both been demonstrated with 
photorefractive screening solitons 

We look for solitary waves of Eqs. (Q) in the form 
u = U(r) exp (i/3 u z), w = W(r) exp (i/3 w z), where the 
envelopes U and W satisfy the equations 
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Here r = ^ x 2 + y 2 is the radial coordinate, and (3 U 
and (3 W are nonlinearity-induced shifts of the propaga- 
tion constants. System (||) has two families of soli- 
ton solutions: U = G u {f3 u ,r), W = and U = 0, 
W = G w ((3 w ,r), which can be found numerically by solv- 
ing the equation G" a + G'Jr - f3 a G a - G a /(1 + G 2 ) = 0, 
where a = {u,w}. These solutions can be characterised 
by the soliton powers P(/3 a ) = 2tt G 2 (f3 aj r) rdr. 

In addition to the one-component solitons, at f3 u = 
f3 v = (3 there exists a family of two-component (vector) 
solitons defined as: U = G((3, r) cos #, W = G(/3, r) sin 6, 
where the variable characterises a power distribution 
between the components. Moving solitons of Eqs. (|l|) 
can be obtained by a well-known gauge transformation. 

To study the soliton collisions analytically, we employ 
a Lagrangian formalism JTo| ? |TTt] ■ Equations (|l|) can be 
obtained from the Lagrangian density: C = (i/2)(u*u z — 
uul) ~ (Kl 2 + M 2 ) + (i/2)(w*w z - ww* z ) - (K| 2 + 
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|i^| 2 ) — ln(l + \u\ 2 + \w\ 2 ). Now we consider the inter- 
action between two spatial solitons (ui,w\) and (u2, ^2), 
taking u = u\ + ^2, w = w\ + W2 and introducing the fol- 
lowing free parameters (j — 1, 2): the positions of soliton 
centers (xj,yj), and the common and relative phases of 
the soliton components uj and Vj , which we denote by 
and i/jj, respectively. Following (ll), we assume that the 
soliton parameters vary slowly in z, and integrate the La- 
grangian density over x and y. After this averaging pro- 
cedure, we reduce the number of equations by using the 
conservation of the angular momentum, Ai = sVbP/4, 
where s is the impact parameter (defined as the minimum 
distance between the trajectories of non-interacting soli- 
tons), and Vo = dRo/dz is the initial value of the soliton 
relative velocity (see Ref. ||). 

The averaged Lagrangian can be presented as L = 
L1+L2 — /7i n t, where the first two terms are the individual 
contributions of the vector solitons, and the third term 
corresponds to an effective interaction potential given by 



tfint = M R s 2 V 2 /(2R 2 ) - U incoh (R) - U coh (R)x 
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where M R = P/2, 0± = 2 ± 6 lt 4>± = <p 2 ± <Ai, 
tp± = ip 2 ± ipi, and R = ^ (x 2 - Xi) 2 + (y 2 - yi) 2 
is the relative distance between the interacting soli- 
tons. The functions Ui nco h and U co h are expressed 
in terms of the soliton overlap integrals, C/ co h = 
2 / !T 00 [G\G 2 /{l + G\) + GlGt/il + G 2 )]dx dy, U incoh = 
I CJG\G 2 /{1 + G 2 ) + G\G\I(\ + G 2 )]dx dy. The first 
term in Eq. ^ describes a centrifugal force (which 
is always repulsive), the second - incoherent attrac- 
tion, and the third - coherent interaction. When R is 
large enough, the soliton interactio n is det ermined by 
the tail asymptotics G(r) ~ exp [— v(lT 0)r]/y/r, that 
yields, U coh (R) ~ exp [- ^(1 + 0)R]/VR, U incoh (R) ~ 
exp [—2^/ (1 + (3)R]/R and U co h ^> fTmcoh- For smaller 
R, although U co h > Uincoh, ^incoh is also important. 

The average Lagrangian generates the following equa- 
tions, 
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where the dots stand for derivatives in z, and the ef- 
fective masses are: M<j, = —dP/d/3, and = Mq = 
(cos 2 0- — cos 2 First, we consider a reduced model 

assuming an additional symmetry, 6+ = tt/2. Then, the 
resulting system has stable stationary points. Solving this 
reduced system numerically we observe linear and even 
strong nonlinear oscillations near the stable minima. In 



general, the period of these oscillations in R is different 
from the periods of 0- and ip- . A stable stationary point 
corresponds to a smooth spiraling of the solitons. 
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FIG. 1. Stable soliton spiraling observed in direct model- 
ing of Eqs. (1) for = -0.5, R = V$00, s = 10, Vo = 0.2, 
and 9- = 0+ — tt/2. (a) 3D view; (b) small oscillations of the 
relative distance between solitons; (c) large-amplitude oscil- 
lations for 0- (quasi-periodic power exchange). 
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FIG. 2. Examples of different dynamical regimes of the soli- 
ton interaction obtained by direct modeling of Eqs. (1) for 
P = -0.5, Vo = 0.2, and fl_ = 6> + = tt/2. The initial sep- 
aration is defined as Ro — V400 + s 2 . (a) Weak soliton in- 
teraction (no spiraling); (b) stable spiraling; (c)(d) unstable 
spiraling (decay of spiraling configuration via resonances). 

However, the analysis of the full dynamical system (^) 
brings a surprise: stable stationary points are absent. 
The main reason for this is the negativeness of the ef- 
fective mass in Eqs. (|j), which is a typical desta- 
bilisation mechanism for any coherent soliton interac- 
tion [||. However, numeriacl simulations show that sta- 
ble dynamical spiraling is still possible. To understand 
the physical mechanism of such a dynamical stabilisa- 
tion, we analyse the effective interaction potential (||). 
Although U co h > Z7i nC oh even for small R, large-scale pe- 
riodic quarter-period out-of-phase oscillations in 0_ and 
ip_ can significantly suppress the effective value of the 
U co h term, thus lowering its maximum value by a factor 
of 5 or more. As a result, the incoherent attraction dom- 
inates and solitons become trapped in a spiraling config- 
uration with oscillations near some i^min and large-scale 
quasi- periodic oscillations in both 6- and ^_ (see Fig. 

i). 
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Solving Eqs. ([l]) [and also (|j)] numerically, we con- 
firm the mechanism of the dynamical soliton spiraling. 
In summary, our theory and numerics show that 

(i) Trapping of two beams in a stable spiraling is possi- 
ble for a large range of parameters [examples are shown, 
e.g., in Fig. 1, for s = 10, and Fig. 2(b), for s = 7]. 

(ii) Initially mutually incoherent colliding solitons [i.e., 
T (0) = 7r/2] become partially coherent due to a periodic 
power exchange between their components. Moreover, 
stable spiraling is always accompanied by a large-scale 
periodic power exchange. 

(hi) Initially introduced partial coherence between in- 
teracting solitons (seed mutual coherence) can result in 
repulsion of out-of-phase solitons and fusion of in-phase 
solitons, preventing spiraling. In this sense, modifying the 
initial mutual coherence can easily transform stable spi- 
raling into repulsion ( "escape" ) or fusion. 

(iv) For smaller s and also for some values of s where 
the spiralling and power-exchange frequencies become 
commensurable, the soliton spiraling is not possible [see 
Figs. 2(c,d)]. A series of 'resonance windows', similar to 
those discovered for 2D soliton interactions fl^l , are ob- 
served. For such values of s, oscillations in R are stronger 
and U co h can become dominant (even being effectively 
suppressed), thus leading to a decay of spiraling. 

To verify our theory, we perform a series of experi- 
ments. The experiments are carried out using the pho- 
torefractive screening nonlinear ity |3|||. In essence, the 
photorefractive nonlinearity is anisotropic |l3[| , which 
makes it non-ideal to test our model. However, many 
experimental results suggest that for a large range of pa- 
rameters, the anisotropy is fairly small: isolated 3D soli- 
tons are almost fully circular and planar collisions 
between 3D coherent solitons are almost fully isotropic 
[jn|, except for a special case, e.g., when the collision 
plane is normal to the c-axis of the crystal and for a par- 
ticular initial distance between the solitons Q. In this 
respect, even though the photorefractive nonlinearity is 
not isotropic in 3D, one can still employ it to qualita- 
tively study the predictions of our theory. We therefore 
extrapolate the known analytic results for 2D photore- 
fractive screening solitons |J, which were all confirmed 
experimentally |l7]], to 3D which concurs with Eqs. ([I]). 



Input 



Output 




(b) l 

(c) l 

(d) l 

(e) 



B A 






bo y 





A' 


B' 










A" 








B • 


a m 





134 urn 



FIG. 3. (a)-(i) Collision of two initially mutually incoher- 
ent solitons with different values of the impact parameter. 
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FIG. 4. (a),(b) Stable soliton spiraling of initially inco- 
herent beams. (c),(d) A "seed coherence" beam C, which is tt 
out of phase with B, prevents the spiraling; (e),(f) When the 
beam C is in phase with B, it causes all three beams to fuse. 

The experimental setup is similar to that of Ref. [||. 
Two soliton beams A and B of wavelength 488 nm, with 
power in the order of /iW, and radii of 12/im FWHM are 
launched into an SBN crystal whose electro-optic coeffi- 
cient is 278 pm/V and the length is 6.5 mm. The initial 
y coordinate of B is 9 jam higher than that of A, and B 
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is launched with its initial trajectory inclined (relative to 
that of A) by an angle of 0.01 radians in the x direction 
and 0.0012 radians in the y direction. The intensity ratio 
between the soliton peak and the background illumina- 
tion, \u\l g = \w\l g , is about 5. A field of 4.2 kV/cm is 
applied against the the crystalline c-axis to generate the 
solitons. The impact parameter is adjusted by shifting 
the initial x coordinate of B at the input, while all other 
initial conditions are kept unchanged. When the separa- 
tion in the x coordinate is larger than 31/im as shown in 
Fig. 3(a)-3(c), solitons A and B barely interact [compare 
with Fig. 2(a), which shows a passing of two solitons as 
if they do not interact at all] . As the impact parameter is 
reduced by shifting B closer to A, as shown in Fig. 3(d)- 
3(f), A and B's trajectories are bent due to the attraction 
force between them, and the amount of bending (scatter- 
ing) is dependent on the impact parameter. This mimics 
a classical particle scattering experiment. We distinguish 
A from B and measure the power exchange by monitor- 
ing the output within a time window much shorter than 
the response time of the SBN crystal (1 sec) after A or 
B is blocked. The measured power exchange is smaller 
than 1% in Fig. 3(a)-3(f). 

When we further reduce the separation in the x coor- 
dinate to 9/im [Fig. 3(g)], the two solitons rotate around 
each other [cf. Fig. 1 and Fig. 2(b)]. We find that 60% 
of A and 46% of B at the input go to A 7 (at the output) 
and the rest goes to B' . This power exchange is what 
we have called induced coherence. We also find that a 
small variation in B's initial position or trajectory which 
does not change the rotation angle of beam trajectories 
by much, can cause the fraction of the exchanged power 
to vary considerably [compare with Fig. 1(c)]. In some 
spiraling cases, as low as a 5% level of power exchange has 
been measured at the output of the crystal. In a similar 
spiraling experiment, but with different initial trajecto- 
ries, we find that the power exchange also depends on 
the intensity ratio. That is, the level of saturation of the 
nonlinearity. In that experiment, 17% power exchange is 
measured when solitons are generated with the intensity 
ratio of 12 and only 2% for the intensity ratio of 4. 

We then reduce the x separation further to 4 jiia [Fig. 
3(h)], and find that A and B interact strongly, but the 
spiraling seems to be unstable [compare to the numerical 
result shown in Fig. 2(c,d)]. Finally, when B is launched 
with its initial position beyond A [Fig. 3(i)], they simply 
escape from each other. 

In order to study how the initial partial coherence af- 
fects the soliton interaction, we introduce at the input a 
"seed coherence" beam C which is coherent with B but 
overlaps entirely and copropagates with A. When C is 
added, the intensity of A is reduced to make the total 
intensity (A+C) equal to that of B. The relative phase 
between C and B is adjusted with a tilted piece of glass. 
Before C is launched, we make sure the initial conditions 



of A and B generate a spiraling pair [Figs. 4(a,b)]. When 
C is first adjusted to be out of phase with B [indicated by 
the dark notch between them at the input, Fig. 4(c)], B 
and A+C cannot spiral but just escape from each other, 
as shown in Fig. 4(d), although the power in C is only 
about 28% of A+C. When C is in phase with B, as shown 
in Fig. 4(e) (each intensity of A or C is 50% of B), A, B 
and C fuse into one beam [Fig. 4(f)]. These experimen- 
tal results agree with our theory, emphasizing the fact 
that seed coherence can be used to control the interac- 
tion outcome: spiraling, "escape" , or fusion. 

In conclusion, we have analysed fully 3D interaction 
and spiraling of spatial solitons in an isotropic saturable 
bulk medium. The analysis, numerical simulations, and 
a series of experiments have revealed the important phys- 
ical mechanism of the stable spiraling: a periodic power 
exchange between the interacting beams via induced co- 
herence. Our results and conclusions are expected to hold 
for other types of (even anisotropic and nonsaturable) 
nonlinearity that depends on the total beam power and 
supports stable self-trapped beams in a bulk. 
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